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Abstract: 

We experimentally demonstrate for the first time that a linearly polarized beam is fo- 
cussed to an asymmetric spot when using a high-numerical aperture focussing system. 
This asymmetry was predicted by Richards and Wolf [Proc. R. Soc. London A 253, 358 
(1959)] and can only be measured when a polarization insensitive sensor is placed in the 
focal region. We used a specially modified photodiode in a knife edge type set up to 
obtain highly resolved images of the total electric energy density distribution at the fo- 
cus. The results are in good agreement with the predictions of a vectorial focussing theory. 

The ability to control and focus light is at the heart of the interdisciplinary field 
of optics and for many applications an exact knowledge of the structure of the focal 
field is required. Many areas in optical sciences make use of a tightly focussed light 
beam such as confocal microscopy ^ and optical data storage . A highly concentrated 
and well matched field is also a necessary requirement for coupling to small quantum 
systems jS] and applying light forces to microscopic particles |3]. In the regime of strong 
focussing the widely used scalar theories are inadequate to describe the focal field. A 
vectorial focussing theory is required instead. In the concrete example of a beam which 
is linearly polarized along the y-axis, the rays propagating in the xz- and in the yz-plane, 
respectively, contribute differently to the focal field (see figures Ufa), (b)). This has two 
closely related consequences. The effective numerical aperture in the yz-plane is reduced 
and the focal spot is elongated in the direction of polarization of the input beam because 
the rays that propagate in this plane do not add up perfectly at the focus. Associated with 
this elongation, field components in directions orthogonal to the direction of polarization 
of the input beam arise. These effects are neglected in scalar descriptions 0, which 
therefore do not predict the correct spot shape and size. 

All aspects of the field distribution of a tightly focussed linearly polarized light beam have 
been studied in theory taking into account the vector properties of the field [HI 13 El IH]- The 
theoretical calculations predict pronounced deviations from the results of the scalar theory. 
Nevertheless, a major part of the intensity^ is still predicted to be contained in a field 

"'^In the following, intensity always refers to electric energy density, which is the part of the field energy 
that couples to standard photodetectors and photosensitive materials. 
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Figure 1: Plane wave model for focussing a beam which is linearly polarized parallel to the 
y-axis. Cross sections for two planes which contain the optical ) perpendicular and 

(b) parallel to the direction of polarization, show that the contribution to the total field of 
the rays propagating under a large angle to the optical axis is weaker in (b). Calculated 
intensity patterns for the three orthogonal field components (c-e) and total intensity 
distribution (f). The direction of polarization is indicated by arrows. (Calculations for 
NA=0.9, Gaussian input beam with a diameter (1/e^) of 94% of the diameter of the 
entrance pupil of the objective. 



component with an almost circular intensity distribution, which is polarized parallel to the 
direction of polarization of the incident beam (see figure ^(c)). The remaining intensity 
appears in the two other mutually orthogonal field components. The orthogonal transverse 
component is calculated to show a four leaf clover shaped distribution (see figure[T](d)) and 
does not contribute substantially to the shape of the focal field due to its minute intensity. 
On the other hand, the intensity contained in a longitudinal z-component is calculated to 
contribute 16% to the total intensity (for a numerical aperture of 0.9 and for homogeneous 
illumination). The calculated pattern of the intensity distribution of this component 
shows maxima away from the optical axis in the direction of the incident polarization 
(see figure H(e)). This causes the asymmetry of the total field (see figure ^(f)). When 
annular apertures are used for focussing, the relative strength of this longitudinal field 
and therefore the asymmetry is predicted to be enhanced 

On the experimental side several studies of the focal intensity distribution have been 
reported. Using a CCD-camera Karman et al. ^1] achieved a large dynamic range of 
5 orders of magnitude for the intensity. They focussed with a low numerical aperture. 
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NA=0.01, for which polarization effects are too small to be detectable. The results were 
in excellent agreement with the predictions of the scalar theory. For the observation of 
polarization effects it is not only important to focus the light with a high numerical aper- 
ture but also to detect light with an equally high effective numerical aperture. Firester et 
al. moved a knife edge through a strongly focussed beam (with numerical apertures 
up to 0.95). They observed only a small solid angle of the transmitted light because the 
detector was not in direct contact with the knife edge. Therefore they did not observe 
any asymmetry. Instead of the knife edge, one can also move a Ronchi ruling through the 
beam but in this case there is an additional requirement: one has to make a priori as- 
sumptions about the shape of the intensity distribution for the evaluation [TH]. Collecting 
the light with a fibre tip typically provides the necessary high numerical detection aper- 
ture Plj. However, the metal coating around the tip introduces a significant polarization 
dependence [THj. Another option is to scan small sub- wavelength beads through the focus 
and observe the scattered and reflected light jTHllIIj. This technique, however, does also 
not readily disclose the focal intensity distribution in an unambiguous way. The detected 
signal depends on the emission characteristics of the bead [18j which in turn depends 
on the direction of the exciting field. This direction varies largely in the case of high 
numerical aperture focussing. If the detection solid angle is not exactly equal to 2n or 47r, 
the overall detection efficiency will depend on the direction of the field exciting the bead. 
Therefore, one again has to make a priori assumptions about the measured distribution 
of the vectorial electric field E. Instead of isotropic beads, which do not have a preferred 
axis, one can also use individual oriented molecules the excitation of which depends sen- 
sitively on the direction of polarization of the incoming light. This property was used 
to determine the orientation of molecules on a surface Single molecules oriented in 
three different mutually orthogonal directions could be used to scan the intensity distri- 
bution of all three polarization directions separately. However, a quantitative comparison 
of the three patterns remains difficult as the correct relative scaling for the composition 
of the total intensity pattern requires the exact knowledge of the polarization-dependent 
sensitivity of the detection system. An independent measurement which relates the po- 
sition integrated intensities of different polarization components has been obtained only 
for the two transverse components |2(]J. To summarize the state of the art, individual 
aspects of the effects of strong focussing have been measured in several experiments, but 
none of them has measured the asymmetry of the polarization integrated total intensity 
distribution in the focal spot. 

We present here the first direct experimental measurement of the total intensity dis- 
tribution summing all three polarization components. We verify the elongated shape of 
the total intensity distribution due to the presence of all three mutually orthogonal field 
components. The results of our measurements are in agreement with the results of the 
vectorial focussing theory [HIIZI- 

To measure the intensity distribution in the focus we use an experimental set up which 
is based on the knife-edge method. In this set up the transmitted light is detected in the 
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immediate vicinity of the knife edge. Thus the effects of diffraction and propagation of 
the beam are minimized. Our sensor is a specially designed p-i-n photodiode that is 
moved transversely in the focal plane using a piezo driven table operating in closed loop 
mode with a step width of lOnm to ensure high lateral resolution |2II. Parts of the light 
sensitive area are covered by a sharp edged opaque pad that plays the role of the knife 
edge. The photocurrent of the diode is proportional to the electric field intensity (oc \E\'^) 
that is incident on the light sensitive area, which depends on the fraction of the beam 
blocked by the opaque structure. A more detailed description of the photodiode is given 
in the appendix. The data were taken with a laser intensity noise <0.6% and the signal 
noise was reduced to <0.3% by normalizing to the signal of the monitor diode. At each 
position of the edge the photocurrent is averaged for up to 100ms which further reduces 
the signal noise to <0.1%. For every scan of the edge the signal is corrected for offsets 
when the beam is completely blocked by the opaque structure and is normalized to one 
when the whole beam hits the light sensitive area. The differentiation of the signal with 
respect to the detector position yields the line integrals of the intensity distribution in one 
transverse direction, parallel to the edge. This is equivalent to the projection of the 2D 
intensity distribution on to a line perpendicular to the edge. The curves are differentiated 
with the Savitzky-Golay algorithm |22] using five adjacent points. 



Figure 2: Experimental set up to measure the focal intensity distribution. MO: Leica 
plan-apo lOOx, 0.9, OD: optical diode, HWP: half-wave plate, TL: telescope lenses, PH: 
pinhole, FL: focussing lens, NCFPI: non-confocal Fabry-Perot interferometer, CL: col- 
limating lens, AS: aperture stop, M: four mirrors for polarization insensitive deflection 
(details not shown), MD: monitor diode, PD: photo diode partially covered with gold-zinc 
alloy 

To provide a linearly polarized and well defined input field, a helium neon laser beam 
(A=632.8nm) is sent through a non-confocal Fabry-Perot interferometer operating as a 
mode cleaner. The transmitted beam is coUimated to a TEMoo-mode with a beam radius 
of ujG=1.7mm (see figure EI). Using a half-wave plate, the direction of polarization was 
rotated with respect to the edge. There was no significant beam displacement due to the 
rotation of the half-wave plate, documented by the fact that the maximum displacement 
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of the peak positions in the sensor plane is below 20nm. Measurements for 18 different 
orientations of the polarization relative to the edge were performed to reconstruct the 
intensity distribution in the focal plane using the Radon back transformation [2H]. For 
measurements with an annular aperture, a stop is placed directly in front of the microscope 
objective to block the inner part of the input beam. The stop is a glass substrate coated 
with an opaque disc (dj=2.0mm). The diameter of the entrance pupil of the microscope 
objective is 3.6mm. 

In figure E^a) and (b) we present distributions of the focal intensity for a linearly 
polarized beam that was focussed without and with annular apertures respectively. The 
distributions were reconstructed from the experimental data and both show the elongation 
of the focus. It is clearly visible that the elongation of the focus increases when an annular 
aperture is inserted. Figures El (a) and (b) show cross sections through the reconstructed 
distributions. The focal spot size is reduced from 0.40A^ (theory: 0.37A^) without annular 
aperture to 0.35A^ (theory: 0.28A^) with annular aperture. This is mainly due to a reduc- 
tion of the width of the profile in the direction perpendicular to the initial polarization. In 
this direction the increase of the sidelobes due to the annular aperture is most pronounced. 
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Figure 3: (a) Reconstructed intensity distribution in the focal plane for a linearly polarized 
input field without annular apperture. The contour line at half the maximum value (thick 
black line) is shown rotated by 90° (dashed line). The spot is elongated in the direction 
of polarization (indicated by the arrow). Using an annular aperture (b) the elongation 
increases. 



The asymmetry is even more obvious if one goes back to the input data for the tomographic 
reconstruction i.e. the initial line integrals of the intensity distribution (see figures |3] 
(c)and (d)). Parallel to the direction of polarization the distribution is broader than in 
the orthogonal direction. 

It is remarkable that the measured and the theoretical profiles are in good agreement 
even in the shapes of the slopes. The measured height ratios of these projections of 0.83 
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(theory: 0.79) without and 0.73 (theory: 0.74) with annular aperture are close to the 
theoretical values. 




p / |jm p / |jm 

Figure 4: Cross sections through the reconstructed intensity distributions (a,b) along a line 
parallel (solid circles) and perpendicular (open squares) to the direction of polarization. 
Values for the measured FWHM are provided together with theoretical values in brackets. 
(c,d) Measured projections of the intensity distribution on to a line parallel (solid circles) 
and perpendicular (open squares) to the direction of polarization. (a,c) without and (b,d) 
with annular aperture. Solid lines: theoretical curves 

Also note, that we carefully checked that the elongation of the focus is not an artefact 
caused by any residual polarization dependent interaction with the edge. The experi- 
mental evidence is that the elongation completely vanishes when the effective numerical 
aperture of the microscope objective (Leica, Plan Apo, lOOx, 0.9) is reduced stepwise to 
0.27 by restricting the diameter do of the input beam by using an aperture stop (figureEI). 

The observed deviation between theory and experiment, 5% for the measured data and 
up to 20% for the deduced focal area, may be attributed to residual imperfections of the 
set-up each causing deviations at the percentage level. Firstly, the Debye-approximation 
used to obtain the theoretical curves neglects diffraction at the aperture rim. Second, 
the objective is man made; in addition to radially dependent Fresnel-losses, there may 
be stress-induced birefringence and deviations from the perfect surface geometry. Third, 
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Figure 5: Peak height of the projected intensity distributions vs. effective numerical aper- 
ture. Sohd squares: projections on to a hne perpendicular to the direction of polarization, 
open circles: projections on to a line parallel to the direction of polarization along with 
theoretical curves (solid lines). An input field with a homogeneous intensity distribu- 
tion was used. Note that the polarization anisotropy appears only for high numerical 
apertures. 



similar arguments hold for the overall beam geometry and the detector. The positioning 
precision e.g. of the annular aperture with respect to the objective is about 50/im. Fourth, 
the Radon back transformation requires perfectly centered data. A centering error of 10~^ 
causes already changes of about 10^^ for the spot area. Fifth, the differentiation of the 
signal is sensitive to beam instabilities. The mechanical stability of the set up has been 
improved until the position fluctuation of the focal spot was below Inm. We estimate 
that the effect of all these imperfections of the set up and of the theory can well explain 
the observed discrepancy. 

Line integrals of the intensity distribution were measured also for different defocus 
positions (see figure E)). This yields a cross section of the propagating beam containing 
the optical axis. The observed asymmetry with respect to the focal plane is caused by 
minute deviations of the wavefront |21 introduced by the microscope objective which are 
at the limit of a state-of-the-art production process. For the comparison of the experi- 
mental results with numerical calculations we considered spherical aberrations and found 
the best agreement for a wavefront with an rms error of 0.00 75 A, which is in excellent 
agreement with the manufacturer's specifications [22] conveyed to us after the measure- 
ment. Note, that the projections in figure 6 are distinctly different for the two planes 
parallel and perpendicular to the initial polarization. Also note, that the effect of these 
aberrations was observed in the defocus series. According to theory it does not influence 
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the focal size and shape in a measurable way. 
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Figure 6: Projections of the intensity distribution for a linearly polarized input beam. 
Measured (a,c) and calculated (b,d) projections on to planes which contain the optical 
axis which are parallel (a,b) and perpendicular (c,d) to the direction of linear polarization. 
A spherical aberration (0.0075A rms) of the wavefront was assumed in the calculations. 



In summary, we report the first reconstruction of the elongated focus caused by polar- 
ization effects. The set up is also sensitive enough to detect quantitatively the effect of 
the minute spherical aberrations introduced by the focussing microscope objective. The 
longitudinal field component showing up in the focal plane may help to improve the spec- 
troscopic investigation of the 3D orientation of a single molecule near a surface [26i ITTj 
and evidence for the existence of the longitudinal field component has already been re- 
ported jini. The good agreement between the measured and calculated intensity profiles 
we report here allows one also to quantitatively infer the appearance of a longitudinal 
field component. 
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Appendix: 

The photodiode used in our experiment is a p-i-n AlGaAs diode that consists of a GaAs 
substrate, an n-doped layer of AlGaAs, a layer of intrinsic GaAs, a p-doped layer of 
AlGaAs and a p-doped GaAs cover layer (seen from the substrate to the surface). The 
cover layer lOnm thick is necessary to protect the AlGaAs from corrosion. During the 
fabrication process, care must be taken not to damage this layer. For a diode that 
showed scratches that were up to 30nm deep as measured by an AFM, the photocurrent 
showed unacceptably large fluctuations when the focussed spot was scanned over the 
active surface. The diode used in the experiment showed a surface roughness of below 
3nm rms and the corresponding fluctuations were below 0.1% (peak-to- valley). The active 
surface of the device is 200/im x 260^m wide and is surrounded by rectangular contact 
pads to tap the photocurrent which each have a size of 200/xm x 120//m. 
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Figure 7: Schematic distribution of the photo current (b) and its derivative (a) with 
respect to the position of the focal spot relative to a grating structure (c) placed on the 
photodiode. The distance between the extrema before and behind one bar of the grating 
in (a) depends on the direction of polarization of the incident beam relative to the grating, 
d^ denotes this distance for a polarization perpendicular (solid line) to the bars and dp 
the distance for a polarization parallel (o) to the bars. 



The opaque material that forms the knife edge plays an important role. To allow for a 
thin layer the opaque structure should show a small transmission and therefore preferably 
be made of metal. However, for a metallic edge material, one would typically expect 
a polarization dependent interaction with the light field. For several edge materials we 
observed that the position of the strongest increase of the photocurrent depends on the 
direction of the polarization of the input field with respect to the edge when scanning the 
edge relative to the focus. To characterize this polarization dependent shift we determined 
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the "duty-cycle" of the signal measured for different grating structures on the active 
surface. Each grating was measured twice once with the polarization of the input field 
being parallel and perpendicular to the grating. The difference in the " duty cycle" for the 
two directions of polarization of the input field allows one to infer the relative lateral shift 
of the two transverse field components in the focal plane (see figure Hj). This procedure 
has the advantage that it does not require any absolute position measurements. 

The shift was determined for different edge materials, gold and different gold-zinc 
alloys. It turned out that this shift depends on the zinc-concentration of the edge material 
(see figure IHl). 



Figure 8: Measured difference (d^— dp) of the interval between the extrema in the line 
integrated projections. This difference is caused by the polarization dependent interaction 
with the edge as a function of the zinc concentration. A least square fit (dashed line) was 
used to find the optimum zinc concentration (about 20%). 

We found no significant relative shift for a zinc concentration of about 20%. Therefore, 
all measurements presented in the main part of the paper were performed using the edge 
material with the optimized concentration and a thickness of about < 200nm, resulting 
in a transmission below 1%. For the chosen alloy a good agreement with the results pre- 
dicted by theory was observed. All other materials which were tested so far were excluded 
as well behaved edge materials. It is worth noting that for the well behaved edge that 
caused a minimal differential shift between the x and y polarization component (see figure 
8) we also did not observe any significant shift for the z-component. Here x and y refer 
to the two transverse components and z to the longitudinal polarization component, x 
is parallel to the edge whereas y and z are perpendicular to the edge. The mechanism 
that for some materials leads to a displacement of individual field components is not yet 
completely clarified. A possible mechanism under study is the excitation of surface plas- 
mons. However, we want to stress that the above mentioned test gives phenomenological 
evidence that the diode with the well behaved edge material is suitable to perform the 
described measurements of the asymmetry of the focus. 
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